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ABSTRACT 

A microscopic analysis of the viscous energy gain of energetic particles in 
(gradual) non-relativistic shear flows is presented. We extend previous work and 
derive the Fokker-Planck coefficients for the average rate of momentum change 
and dispersion in the general case of a momentum-dependent scattering time 
r(p) oc p a with a > 0. We show that in contrast to diffusive shock acceleration the 
characteristic shear acceleration timescale depends inversely on the particle mean 
free path which makes the mechanism particularly attractive for high energy seed 
particles. Based on an analysis of the associated Fokker-Planck equation we show 
that above the injection momentum p power-law differential particle number 
density spectra n(p) oc are generated for a > if radiative energy losses 

are negligible. We discuss the modifications introduced by synchrotron losses 
and determine the contribution of the accelerated particles to the viscosity of 
the background flow. Possible implications for the plasma composition in mildly 
relativistic extragalactic jet sources (WATs) are addressed. 

Subject headings: acceleration of particles — galaxies: jets 

1. Introduction 

Shear flows are a natural outcome of the density and velocity gradients in extreme as- 
trophysical environments, and in the case of active galactic nuclei (AGNs), for example, 
observationally well-established (cf. Laing & Bridle 2002; Laing et al. 2006; Rieger & Duffy 
2004). The acceleration of energetic particles occurring in such flows can represent an effi- 
cient mechanism for converting a significant part of the bulk kinetic energy of the flow into 
nonthermal particles and radiation as has been successfully shown in a number of contribu- 
tions (e.g., Ostrowski 1990, 1998; Stawarz & Ostrowski 2002; Rieger & Duffy 2004, 2005c). 
Early theoretical progress in the field has been achieved based on a detailed analysis of the 
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Boltzmann transport equation (Berezhko & Krymskii 1981; Earl et al. 1988; Webb 1989). 
Somewhat similar to the microscopic picture for Fermi acceleration, shear acceleration is 
essentially based on the fact that particles can gain energy by scattering off (small-scale) 
magnetic field inhomogeneities moving with different local velocities due to being embed- 
ded in a collisionless shear flow with a locally changing velocity profile (see Rieger & Duffy 
2005b, for a recent review). In a scattering event particles are assumed to be randomized 
in direction with their energies being conserved in the local (comoving) flow frame. As the 
momentum of a particle travelling across a velocity shear changes with respect to the local 
flow frame, scattering leads to a net increase in momentum with time for an isotropic particle 
distribution. In contrast to 2nd order Fermi acceleration, however, not the random but the 
systematic velocity component of the scattering centres is assumed to play the important 
role. 

In the present paper we present a microscopic analysis of shear acceleration where the under- 
lying physical picture becomes most transparent. Sect. 2 gives the derivation of the Fokker- 
Planck coefficients for a simple, non-relativistic (longitudinal) gradual shear flow. Sect. 3 
analyzes the corresponding Fokker-Planck transport equation providing time-dependent and 
steady-state solutions for specific cases. The contribution of the accelerated particle to the 
viscosity of the background flow is determined in Sect. 4, while implications for mildly rela- 
tivistic jet sources are addressed in Sect. 5. 

2. Derivation of Fokker Planck coefficients 

For a gradual non-relativistic shear flow, Jokipii & Morfill (1990) have calculated the 
Fokker-Planck coefficients using a microscopic treatment restricted to a mean scattering 
time r c = const, independent of momentum (cf. also Rieger & Duffy 2005a). However, as 
under realistic astrophysical conditions the scattering off magnetic turbulence structures is 
momentum-dependent (e.g., gyro-, Kolmogorov- or Kraichnan-type) , we extend this analy- 
sis to the more general case where the local scattering time r c is a power-law function of 
momentum, i.e., r c {p) = r p a , a > 0. 

We consider a simple non-relativistic two-dimensional continuous shear flow with velocity 
profile given by u = u z (x) e z . Similarly as in Jokipii & Morfill (1990), we choose a spherical 
coordinate system where 9 denotes the angle between the a>axis and the velocity vector 
v = (v x ,v y ,v z ) of the particle, and is the azimuthal angle such that tan0 = v y /v x . Let m 
be the relativistic mass and pi the initial momentum of a particle relative to its local flow 
frame and £>2 the corresponding momentum after the next scattering event. As the particle 
travels across the shear, its momentum (and thus its mean free path) changes with respect 
to the local flow frame. Denoting by f =< r > the mean scattering time averaged over 
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momentum magnitude, one obtains 

r = Tc+ o^ Ap = r c l + -a , 1 

in the limit Ap/pi = (p 2 — Pi)/pi <C 1 where collisions are assumed to produce only a 
small change in the momentum of particles, with large changes occurring only as a conse- 
quence of many small changes. Within the time f a particle travels a distance 5x = v x f = 
(pi/m) cos 9 f in the x-direction across the shear and the flow velocity changes by an amount 
Su = duel, where 

« = (£)« = ft ,( 1 + i a £), (2) 

and Su := (du z /dx) (pi/m) cos9t c . Contenting ourselves with a Galilean transformation for 
the non-relativistic flow speeds involved, the particle's momentum relative to the flow will 
thus have changed to p 2 = pi + m5u, i.e., one obtains 

p\ = Pi + 2 m 5u pi tZ + m 2 (8u) 2 

= p? 1 + 2 sin cos 0+ 5—1. (3) 

V pi pi J 

As the next scattering event is assumed to preserve the magnitude of the particle momentum 
relative to the local flow frame, the particle magnitude will have this value in the local frame. 
To second order in 5u, one thus finds 

1 H Su sin 9 cos H 5- Su 2 

Pi Pi 

_I HL Su 2 (1 - a) sin 2 9 cos 2 A . (4) 

2 Pi / 

By using spherical coordinates and averaging over all momentum directions for an isotropic 
particle distribution, the Fokker-Planck coefficients (Chandrasekhar 1943) describing the 
average rate of momentum change and the average rate of momentum dispersion (associated 
with a broadening of the distribution) can be determined as follows 

Ap\ 2<p 2 - Pl > (4 + a) (du z \ 2 

p hr 7 > ( 5 ) 



At / t c 15 V dx , 

/(Ap) 2 \ 2<(p 2 - Pl ) 2 > 2 2 (du^ 2 



\ At I ' t c 15 F \dx 



^P 2 ~^ r c , (6) 



where the index 1 has been dropped on the right hand side, revealing that a net increase 
in momentum proportional to the square of the flow velocity gradient occurs with time. 



-4- 



Obviously, the stronger the shear and the larger the particle mean free path, the higher 
the possible impact of scattering and thus the higher the rate of momentum change. For a 
constant scattering time, i.e., for a — 0, these expressions reduce to those found by Jokipii 
& Morrill (1990), apart from the factor of two that takes the time average properly into 
account. The Fokker-Planck coefficients are related by the equation 

/ Ap\ _ j_ d_ ( 2 /(Apf\\ _rfl, M 

\ At / ~ 2p 2 dp \ P \ At I)' p*d P [P Tc) ■ [ ' 

where V on the right hand side denotes the shear flow coefficient (cf. Earl et al. 1988) which, 
for the flow profile chosen, is given by 

1 / dui du A 2 2 duiduj 1 f du z \ 2 

on + ~aZT ~ 77 o_ o - T7 "aT ■ W 



30 (9xj y 45 cte; ctej 15 \ 9x _ 

Eqs. (5) and (7) indicate that for the adopted scaling r c oc p a acceleration occurs as long as 
a > —4. Comparison with previous approaches shows that the average rate of momentum 
increase < Ap/At >= (4 + a)Tpr c , as given by eq. (5), agrees with the non-relativistic 
limit of the full shear acceleration coefficient derived in Rieger & Duffy (2004) and the re- 
sults found by Berezhko & Krymskii (1981) and Earl et al. (1988). 

Eq. (5) implies a characteristic acceleration timescale for particle acceleration in non-relativistic 
gradual shear flows 

t = E = °- (9) 

acc < Ap/At > (4 + a)TA' [ ) 

which is inversely proportional to the particle mean free path A ~ cr c , in remarkable contrast 
to diffusive shock acceleration where t acc oc A (e.g., Kirk & Dendy 2001). As the mean 
free path of a particle increases with energy it follows that higher energy particles will be 
accelerated more efficiently than lower energy particles. This characteristic behaviour makes 
shear acceleration particularly attractive for the acceleration of high energy seed particles 
and the production of UHE cosmic rays (e.g., Rieger & Duffy 2005c). In the case of a particle 
mean free path proportional to the gyroradius r g , for example, the acceleration timescale 
scales with the particle Lorentz factor in the same way as the timescales for synchrotron and 
inverse Compton losses, respectively. 



3. Particle transport and power-law formation 

The propagation of energetic charged particles in non-relativistic shear flows can be cast 
into a Fokker-Planck-type momentum diffusion equation (cf. Melrose 1980). Taking syn- 
chrotron losses into account the isotropic phase space distribution function f(p,t) (averaged 
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over all momentum directions) then satisfies 



dt p 2 dp \ \ \ At 

+ ^if^) fM ) +liM ' (10) 

where Q(p,t) denotes the source term, < p > s is the synchrotron loss term given by 

<p> s =-f3 s p\ (11) 

with f3 s = 45 2 e 4 /(9w 4 c 6 ) when expressed in Gaussian units, and where the coefficients are 
related by eq. (7). In the absence of losses and injection, eq. (10) reduces to a momentum 
space diffusion equation 

% = (12) 



dt p 2 dp \ dp / 

where D = Tp 2+a T is the momentum space diffusion coefficient. Eq. (7) is known as the 
principle of detailed balance (e.g., Blandford & Eichler 1987) and is a condition which must 
be satisfied for the Fokker-Planck equation to reduce to the case of pure diffusion. 
It is interesting to note that second order Fermi acceleration as a result of scattering from 
forward and reverse Alfven waves is also described by momentum space diffusion (Skilling 
1975). It is also possible to analyse such a process on the microscopic level, deriving the 
Fokker Planck coefficients and demonstrating the principle of detailed balance, see Duffy & 
Blundell (2005) for details. 

In the following eq. (10) will be used to analyze the the evolution of the particle distri- 
bution function for some particular cases. 



3.1. Time-dependent solutions for impulsive source 

Consider first the case where radiative synchrotron losses are negligible and particles are 
injected monoenergetically at time t = with momentum p , i.e., Q(p,t) = QS(p — p ) 5(t). 
Eq. (10) then reduces to 

»fM - 1 a f rA W)\ te((W _ Al , (13) 



dt p 2 dp \ dp 



The same equation is obtained if one uses the generalized Parker transport equation (Earl 
et al. 1988, eq. 6) for the shear flow profile under consideration. As shown in the appendix 
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the mathematical solution of eq. (13) for a ^ takes the form, see eq. (A6), 

fM = <*fl (-Y 3+ ° ,/2 exp Wl ( ' pip) , (14) 

where is the modified Bessel function of the first kind (cf. Abramowitz & Stegun 1972). 
The corresponding differential particle number density p 2 f is illustrated in Fig. (1) using a 
momentum dependence of a = 1. Obviously, the position of the peak shifts to the left 
and decreases with time (thus ensuring particle number conservation) while the distribution 
becomes broader and the relative strength of its tail increases (as a consequence of dispersion 
and acceleration). For small z — > the modified Bessel function is known to approach 
I v {z) ~ (z/2) v /Y{v + 1) for v ^ —1, —2, .. with T(y + 1) denoting the Gamma function. It 
thus follows that for a > and p ~ p the phase space distribution approaches a power-law 
/ oc p~( 3+a ) on a characteristic timescale t c = l/(a 2 rr Po) (cf. Berezhko 1982). For p > p 
the time-integrated solution of eq. (13) will thus show the same power- law behaviour, cf. 
eq. (19). 

Note that the Fokker Planck description employed is only valid at late and not at very 
early stages since the resultant transport equation is essentially non-causal, i.e., eq. (10) 
is a parabolic partial differential equation. Its solutions thus share with those of the heat 
equation the mathematical properties of an infinite speed of propagation (i.e., the distribution 
is non-zero everywhere for t > 0), smoothing of singularities and decrease and broadening 
with time. In general, the problem of infinite propagation may be overcome by using a 
modified diffusion equation of the (hyperbolic) telegrapher's type (e.g., Cattaneo 1948; Morse 
& Feshbach 1953). Formally, eq. (10) relies on the diffusion approximation and therefore only 
provides an adequate description of the particle transport for times much larger than the 
mean scattering time, i.e., as long as r(p) t which constrains the maximum momentum 
range for a given time t over which the solution eq. (14) may be considered appropriate. 
Complementary, from a physical point of view particle energization is expected to occur 
only on a characteristic timescale determined by the systematic acceleration in eq. (9) (e.g., 
Ball et al. 1992), which constrains the characteristic maximum momentum p(t) achievable 
within a given time interval. 

3.2. Time-integrated (steady-state) solutions 

Defining the time-integrated phase space distribution function 

POO 

F(p)= / f(p,t)dt (15) 
Jo 



-7- 



and integrating the Fokker Planck equation (10) over time assuming a monoenergetic source 
term Q(p, t) = Q 5(p - p ) S(t) gives 

Wp ( rp4+v °^r) + h (fo* F ®) = - Q v« 5{ p - p ^ ■ ^ 

Note that the time-integrated function F(p) is proportional to the steady state solution of 
eq. (10) for continuous injection since eq. (16) coincides with eq. (10) if df/dt = 0. The 
general solution of eq. (16) is of the form 

On 2 , s f p p*( p ') 

F(p) = ^ [d - H{p - p )\ e -xW / dp'— + c 2 e~^ (17) 

J p P 

where H{p) is the Heaviside step function, Ci,c 2 are constants to be determined by the 
boundary conditions, and 

. p -( a -i) for « ^ 1 



*&>) = < '- I,r ™k lnp Z'aL'l- <18) 



If synchrotron losses are negligible (f3 s = 0) the time-integrated phase space distribution 
follows a simple power law above p with momentum index — (3 + a) (cf. Rieger & Duffy 
2005a), i.e., one finds 

for a > and boundary conditions F(p — > oo) — > and dF/dp — > for p — > 0. Thus, in the 
case of a mean scattering time scaling with the gyroradius (i.e., Bohm case: a — 1), for exam- 
ple, this corresponds to a differential power law particle number density n{p) oc p 2 F{p) oc p~ 2 
above p an d n (p) oc below p - 

If losses are non-negligible, the solutions become more complex. For a < 1 we end up with 
expressions involving exponential integrals (cf. Abramowitz & Stegun 1972). The charac- 
teristic evolution of the phase space particle distribution function is illustrated in Fig. (2). 
For a < 1 the acceleration efficiency is generally constrained by synchrotron losses, which 
results in a cut-off at p max = ([4 + a] T to//^)^ 1-0 ^ where the acceleration timescale t acc 
(cf. eq. [9]) matches the cooling timescale t syn = l/(j3 s p). If p p max cooling effects are 
negligible and the particle distribution follows a power-law as described in eq. (19). For 
a = 1 the ratio t acc /t syn becomes independent of momentum. Hence if conditions are such 
that (3 s /(Tt ) < 5 (cf. eq. [10]) is satisfied, particle acceleration is no longer constrained by 
synchrotron losses and for (3 S / (T r ) < 4 the resultant particle distribution above p becomes 
F{p) — QPo/ (4 r r — (3 S ) p~ 4 as expected from eq. (19). For a > 1 shear acceleration becomes 
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possible (and is essentially unconstrained by synchrotron losses) when particles are injected 
with momenta above a threshold p m - m = (/3 s /([4 + a] T To)) 1 ^" -1 '. For p ^> p m - m cooling 
effects are completely negligible and the particle distribution approaches a power-law with 
momentum index as given in eq. (19). 

Note that if losses are unimportant a proper physical steady state situation may be achieved 
in instances where particles are continuously injected with momentum po and considered to 
escape above a fixed momentum p max ^> p where the associated particle mean free path 
becomes of order of the size of the system. Formally, this may be done by adding a simple 
momentum-dependent escape term Q(po/Pm&x) 2 8(p—p max ) on the right hand side of eq. (16). 



4. Cosmic-ray viscosity 

As particles gain energy by scattering off inhomogeneities embedded in a background 
flow, efficient shear acceleration essentially draws on the kinetic energy of that flow. In prin- 
ciple, efficient cosmic ray acceleration may thus cause a non-negligible deceleration even for 
(quasi-collimated) large-scale relativistic jets, although a proper assessment of its significance 
in comparison with other mechanisms (e.g., entrainment) will require detailed source-specific 
modelling. In any case, for the purposes of our (non-relativistic) analysis here, the resultant 
dynamical effects on the flow can be modelled by means of an induced viscosity coefficient 
i] s > that describes the associated decrease in flow mechanical energy per unit time, e.g., 

E kin = - Vs J (20) 

for our case of a (non-relativistic) two-dimensional gradual shear flow u z (x) e z (Landau & 
Lifshitz 1982, §16). We may determine this viscosity coefficient using E kin = —E cr , where 
E C t — J fcrdV is the energy gained by the cosmic ray particles (cf. also Earl et al. 1988). 
If particle injection is described by a continuous source term Q5(p — po), the density e cr of 
power gained becomes 

e cr = 4tt J™ p" E(p) dp _ 4 vr Q E(p ) p 2 , (21) 

where E(jp) denotes the relativistic kinetic energy. Note that due to the constraints discussed 
at the end of § 3.1, special care has to be taken if one wishes to evaluate the integral 
in eq. (21), as done in Earl et al. (1988), by replacing the time-derivative of / with the 
Fokker-Planck expression of eq. (13) for the chosen source term. 1 From an astrophysical 



1 Note that as shown above, we cannot choose independently a power-law index for / (e.g., as suggested 
in Earl et al. 1988) once we have chosen a momentum index for the particle mean free path. 



- 9- 



point of view, we are most interested in steady-state type situations where particles are 
injected quasi-continuously with momentum p and, in the absence of significant radiative 
losses, considered to escape above a momentum threshold p max at which A(p max ) becomes 
larger than the width of the acceleration region. The density of power gained then becomes 
e cr ~ AnQplc (jo m axM) - 1) — ^QPoPuva^ c for p max > p > m c 2 . This implies a viscosity 
coefficient 

3 a , . . . 

Vs ^ — KPo) ^oPmax , (22) 

15 

for a > 0, where X(po) — T c (po)c denotes the mean free path for a particle with momentum 
Po, and n = 47r J Pmax p 2 f(p)dp is the number density of cosmic-ray particles in the acceler- 
ation region. 



5. Applications 

Efficient shear acceleration of cosmic-ray particles is likely to occur in a number of 
powerful jet sources including galactic microquasars and extragalactic FR I and FR II sources 
(e.g., Stawarz & Ostrowski 2002; Laing & Bridle 2002; Laing et al. 2006). Whereas a proper 
analysis of powerful AGN-type jets requires a fully relativistic treatment (cf. Rieger & Duffy 
2004), application of the results derived above may allow useful insights in the case of 
moderately relativistic jet sources. As an example, let us consider the possible role of shear 
acceleration in wide angle tailed radio galaxies (WATs), deferring a detailed discussion of 
particle acceleration in microquasar jets to a subsequent paper. WATs are central cluster 
galaxies and appear as "hybrid" sources (Jetha et al. 2006) showing both FR I and FR II 
morphologies. Their inner jets extend tens of kiloparsec, seem to have central speeds in the 
range Vj ~ (0.3 — 0.7) c (provided Doppler hiding effects can be neglected), are apparently 
very well collimated on the kiloparsec scale and exhibit spectra close to S(v) oc z/~°- 5 , with 
evidence pointing to a steeper spectrum sheath likely to be induced by the strong interactions 
between the jet and its environment (e.g., Katz-Stone et al. 1999; Hardcastle 1999; Hardcastle 
et al. 2005; Jetha et al. 2006). Phenomenological studies suggest a velocity transition layer 
width for the large-scale jet of Ar = £ kpc (£ ~ 0.5) and a fiducial jet magnetic field strength 
of B = 10~ 5 &o Gauss (bo ~ 1 for equipartition). As a likely scenario let us consider the case 
where energetic seed particles required for efficient shear acceleration are provided by internal 
shock-type (Fermi I) processes operating in the jet interior and giving rise to a power-law 
proton distribution /<g(p) = UqPi / \i.Tip i ) in the momentum range p\ ~ 2m p c < p < p 2 , 
where p 2 pi may be determined either from the condition of lateral confinement (e.g., 
Pi ~ 10 10 £&o m pC for A ~ Tg) or via the balance of radiative synchrotron losses, and where 
tiq is the number density of accelerated particles. The minimum (Bohm diffusion) timescale 
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for non-relativistic shock acceleration is of the order t acc (p) ~ 6r g c/u 2 where u s is the 
shock speed as measured in the upstream frame (e.g., Rieger & Duffy 2005c). For typical 
parameters u s ~ O.fc, A ~ r g and a simple linear shear profile with V ~ (vj/ Ar) 2 /15, shear 
acceleration (cf. Eq. [9]) will dominate over shock- type processes for protons with Lorentz 
factors above 7 C ~ 10 9 boC,, so that protons with momenta p > p c = ^c^pC may be considered 
as being effectively injected into the shear acceleration mechanism, resulting in a rate of 
injected particles of Q ~ n s /t , where n s ~ n Q (p 1 /p c ) ~ 2 • 10~ 9 n Q , t ~ t acc (p c ) / P csc and 
P csc ~ 4m 2 /c is the escape probability of a particle from the shock (with u 2 measured in 
the shock frame). This yields a viscosity coefficient r/ s ~ n s p max c/(15rt )- Now, the Navier- 
Stokes equation imply that in a viscous flow of density pf and viscosity coefficient r) s , the 
characteristic viscous damping timescale T for the decay of velocity structures of size L is 
of the order T ~ pfL 2 /i] s (e.g., Earl et al. 1988). For L ~ Ar the numbers estimated above 
give a characteristic decay timescale T of the order 

Requiring T to be larger than the timescale t\ set by the apparent stability of the jet flow, 
i.e., U ~ Lj/vj ~ 2 • 10 5 yr for a typical jet length Lj ~ 30 kpc and Vj ~ 0.5 c, gives a 
minimum density ratio contrast of cold matter to energetic protons (Pf /[nQm p }) ~ 10 (0.3/£). 
Therefore, if similar to the supernova remnant or solar wind case atmost a fraction ~ 0.01 
of the overtaken thermal particles is injected into the shock acceleration process (e.g., see 
Trattner & Scholer 1991; Duffy et al. 1995; Baring et al. 1999), significant velocity decay 
effects are unlikely even if the total power of WAT jets should be still dynamically dominated 
by a cold (thermal) proton component. We note that such a case (slight thermal proton 
dominance) appears not impossible and may indeed be consistent with the notion that WATs 
are intermediate stages in a broader framework where FR I appear as pair-dominated sources 
and FR II jets as mainly composed of electrons and protons (e.g., Celotti 1997), although 
circumstantial evidence based on jet bending seems to suggest that WAT jets are rather light 
(Hardcastle et al. 2005; Jetha et al. 2006). 



6. Conclusions 

Turbulent shear flows are widely expected in astrophysical environments. Using a mi- 
croscopic analysis, we have shown that, in the absence of strong synchrotron losses, the 
acceleration of energetic particles occurring in such flows can give rise to power-law dif- 
ferential particle number densities n(p) oc above the injection momentum p for a 
scattering time r c oc p a , a > 0. Dependent on the details of the underlying turbulence spec- 
trum, shear acceleration may thus allow for different power-law indices. As efficient shear 
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acceleration generally requires sufficiently energetic seed particles, this implies an interesting 
corollary: if energetic seed particles are provided by shock-type acceleration processes, the 
take over by shear acceleration may reveal itself by a change of the power-law index above 
the corresponding energy threshold. Perhaps even more interesting is the fact that the char- 
acteristic timescale for particle acceleration in gradual shear flows is inversely proportional 
to the particle mean free path, i.e., t acc oc 1/A. Shear acceleration thus leads to a preferred 
acceleration of particles with higher magnetic rigidity. Indeed, detailed analyses show that 
in realistic astrophysical circumstances efficient shear acceleration works usually quite well 
for protons, but appears restricted for electrons (e.g., Rieger & Duffy 2004, 2005c). As 
shear acceleration essentially draws on the kinetic energy reservoir of the background flow, 
the associated viscous drag force can, depending on the intrinsic plasma characteristics, sig- 
nificantly contribute to a deceleration of the large-scale jet flow. This suggests that shear 
acceleration may have important implications for our understanding of the acceleration of 
cosmic ray particles, the plasma composition and the velocity evolution in astrophysical jet 
sources. 

Financial support through a Marie-Curie and a Cosmogrid Fellowship and discussions 
with John Kirk are gratefully acknowledged. Useful comments by the anonymous referee 
that strengthen the application are appreciated. 

A. Derivation of time-dependent solutions for impulsive sources 

For an impulsive source where Q particles are assumed to be injected with momentum 
p at time t — 0, the shear Fokker-Planck equation reads 



For a^Owe may choose new variables t := Tr t and x := p a , for which the homogeneous 
part of eq. (Al) becomes 



Eq. (A2) can be identified with the Kepinski partial differential equation (Kepinski 1906) 




(Al) 




(A2) 



/ m + ldf ndf 




(A3) 



dx 2 



x dx x di 
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for n = 1/a 2 and m = —(3 + a) /a. For the initial condition f(i — 0, x) — f(x) the solution 
of the Kepinski partial differential equation is known and of the form (Kepinski 1906) 



n 



K a\ m/2 / x + X 
exp I — n 



x 



t 



2n- 



/ xX 



f(X)dX , 



(A4) 



where I v (z) denotes the modified Bessel function of the first kind (cf. Abramowitz & Stegun 
1972). For our initial condition Q5(p — p ) at t = we have 



f(X) = Q\a\p^ a+1) 5(X-X ). 



(A5) 



where A = p a - Using the original set of variables (jp,t), the full solution of eq. (Al) for 
q^O thus becomes 



f(p,t) = 



Qpo 



(a+l) 



Po\ 

|a|rr t \p ) 



(3+a)/2 



exp 



g ° + Po' 
a 2 rrot 



'|l+3/a| 



a 2 Tr p^t 



p_ 

.Po 



-a/2 s 



■ (A6) 



Eq. (A6) agrees with the solution presented in Berezhko (1982) and can be shown to reduce 
to f(p, t) — > Q5(p — p ) in the limit t — > + . In the case of a = — 2 we may use that 



2 \/2nz 



(< 



to obtain 



f(p,t) 



Qpo 



Py/A-KTT t 





(P-Po) 2 ' 




(P + Po) 2 ' 


^exp 


4Tr t 


— exp 


ATT t 



(A7) 



(A8) 



which agrees with the solution derived for this particular case by Earl et al. (1988), see their 
eq. (13). 

For a = we may use the variables z := \np+3t, where t := Tr t, to obtain the characteristic 
diffusion equation d 2 f/dz 2 = df/dt for the homogeneous part of eq. (Al) with fundamental 
solution f(z,t) = exp(— 2 2 /[4t])/(47rt) 1 / 2 . Thus, for a = the solution of eq. (Al) becomes 
(cf. also Kardashev 1962) 



f(p,t) 



Q 



PoVaIFTtoI 



exp — 



In — + 3rr t 
. Po 



/(4rr 



(A9) 
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Fig. 1. — Characteristic evolution of the differential particle number density n(p,t) oc 
p 2 f(p,t) for a = 1 (cf. eq. [14]) as a function of momentum p/po in the case of impul- 
sive injection of particles with po at t = 0. The distribution function is plotted at three 
different times t where t — t' • t c with t c = 1/(TtqPq). The inlay shows the same in double 
logarithm representation and already indicates the formation of a power-law tail n(p, t) oc p~ 2 
for H > 0.3. 
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Fig. 2. — Evolution of the normalized phase space particle distribution function f(p) as a 
function of momentum p/po for different power indices of the mean scattering time r oc p a . 
The thin-dotted lines are drawn to guide the eyes and correspond to power-law distributions 
f(p) p~ 3 ' 5 and p~ 5 , respectively. For a = 0.5 the maximum particle momentum, at which 
acceleration is balanced by losses, has been chosen to be p max = 1000 p , whereas for a = 2 
a minimum momentum j9 min = p /2 has been used. 



